Boundary criticality at the Anderson transition between a metal and 
a quantum spin Hall insulator in two dimensions 
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Static disorder in a noninteracting gas of electrons confined to two dimensions can drive a contin- 
uous quantum (Anderson) transition between a metallic and an insulating state when time-reversal 
symmetry is preserved but spin-rotation symmetry is broken. The critical exponent v that charac- 
terizes the diverging localization length and the bulk multifractal scaling exponents that characterize 
the amplitudes of the critical wave functions at the metal-insulator transition do not depend on the 
topological nature of the insulating state, i.e., whether it is topologically trivial (ordinary insulator) 
or nontrivial (a Z 2 insulator supporting a quantum spin Hall effect). This is not true of the boundary 
multifractal scaling exponents which we show (numerically) to depend on whether the insulating 
state is topologically trivial or not. 

PACS numbers: 73.20.Fz, 71.70.Ej, 73.43.-f, 05.45.Df 



I. INTRODUCTION 

It has long been known that the metallic state of a two- 
dimensional gas of noninteracting electrons is robust to 
sufficiently weak static disorder when time-reversal sym- 
metry (TRS) is preserved but spin-rotation symmetry 
(SRS) is brokcn.El Increasing the disorder strength rela- 
tive to the Fermi energy induces a continuous quantum 
(Anderson) transition to an insulating state. Bulk prop- 
erties of this metal-insulator transition have been well 
characterized numerically when the insulating state is 
topologically trivial.oBu 

It has only been realized in the last few years that there 
are two distinct classes of time-reversal invariant (band) 
insulators: Z 2 topologically trivial and nontrivial insu- 
lators. The Z 2 topological insulator in two dimensions 
supports a helical edge state which is a Kramers pair of 
counter propagating gapless excitations.uu This helical 
edge state is responsible for the quantum spin Hall (QSH) 
effect: an electric field induces a spin accumulation on 
edges transverse to the direction of the electric field 
The QSH effeci..iias been observed in HgTc/(Hg,Cd)Tc 
quantum wells .On Just as edge states in the integer quan- 
tum Hall effect are stable against disorder, the helical 
edge state in a Z 2 topological insulator survives impurity 
scattering as long as Xhe bulk energy gap is open and 
the TRS is preserved. EjE9 This implies that, even in the 
presence of disorder, time-reversal invariant insulators 
can be classified into two distinct classes, Z 2 topologi- 
cal and nontopological (ordinary) insulators, according 
to the presence or absence of a helical edge state. 

It is then natural to ask0 whether the transition be- 
tween the metallic and the QSH insulating state belongs 
to a universality class different than that of the (ordinary) 
two-dimensional symplectic universality class discovered 
in Rcf. § It was shown in Ref. |l2| that the answer is 
negative for the scaling exponent v of the diverging lo- 



calization length upon approaching the Anderson tran- 
sition from the insulating sides. In this paper, we show 
numerically-, that there are boundary multifractal scaling 
exponents^ that are sensitive to the presence or absence 
of a helical edge state on the insulating side of the tran- 
sition. 

We review in Sec. [0] the definition of the two- 
dimensional network model introduced in Rcf. |l2| that en- 
codes the transition between the metallic and QSH insu- 
lating state, as well as the transition between the metal- 
lic and ordinary insulating state in the two-dimensional 
symplectic universality class. The latter transition is 
conventionally studied using the two-dimensional tight- 
binding model introduced in Ref. ^, which is also briefly 
reviewed. The phase diagram for the network model is 
reviewed, and the relevance of boundary conditions to 
the presence or absence of helical edge states is discussed 
in Sec. III. The dependence of the localization length on 



transverse boundary conditions in quasi-one-dimensional 
geometries is discussed in Sec. IV. Boundary multifrac- 



tal scaling exponents in the network model are calculated 
numerically in Sec. 0. Corner multifractal scaling expo- 
nents in the network model arc in vesti gated analytically 
in Sec. 



VI. We conclude with Sec. VII 



II. 



NETWORK MODEL 



Our starting point is a network model introduced 
in Ref. [l2] to capture the Anderson transition between 
the two-dimensional metallic and topological insulating 
states. The network model is constructed by decorating 
an underlying square lattice of sites and single bonds con- 
necting nearest-neighbor sites with the elementary build- 
ing blocks from Fig. [I]. By taking advantage of the bipar- 
tite nature of the square lattice, one colors all sites from 
one sublattice in red and all sites from the complemen- 
tary sublattice in blue. One then replaces any red (blue) 
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FIG. 1: (Color online) Elementary building blocks of the net- 
work model. A square Bravais lattice with nearest-neighbor 
sites connected by bonds underlies the construction of the 
network model. Any red site of one of the two sublattices of 
the square lattice is replaced by a red circle S (a node of type 
S) with the four bonds meeting at the site replaced by four 
pairs of directed links numbered according to the rule shown 
in (a). Any blue site of the complementary sublattice of the 
square lattice is replaced by a blue circle S' (a node of type 
S') with the four bonds meeting at the site replaced by four 
pairs of directed links numbered according to the rule shown 
in (b). Observe that a clockwise rotation by 7r/2 turn (a) into 
(b). The directed links represent incoming or outgoing plane 
waves with a well-defined projection of the spin-1/2 quantum 
number along the quantization axis. Each pair of links re- 
placing a bond represents a Kramers doublet of plane waves. 
Each node S or S' depicts a scatterin g pr ocess represented by 
a 4 x 4 unitary matrix defined in Eq. (2.1) that preserves time- 



reversal symmetry (TRS) but breaks spin-rotation symmetry 
(SRS). 



site with a node S (S') represented graphically by a red 
(blue) open circle. Second, any single bond connecting a 
pair of nearest-neighbor sites of the square lattice is re- 
placed by a pair of directed links of opposite orientation. 
On the links, the spin-1/2 is a good quantum number. A 
link represented by a full line carries the spin-1/2 quan- 
tum number a =f. A link represented by a dashed line 
carries the spin-1/2 quantum number a =[. Third, the 
four pairs of directed links that meet at a node arc la- 
beled according to the rules of Fig. |l|(a) and Fig. |l](b) 
if the node is of type S and S', respectively With the 
conventions of Figs. [j](a) and |l](b) either node defines a 
4x4 scattering matrix S that preserves TRS but breaks 
the SRS and can be represented by 
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cosh X ' 

Q = ia cos 9 sin ip 1 + a l sin 9 cos tp 2 
— a 2 sin 9 sin ip 2 + a 3 cos 9 cos ip 1 . 

Here, the four matrices <t i 2 3 ac t on the spin-1/2 com- 
ponents with (Tq the unit 2x2 matrix and (a 1 , cr 2 , cr 3 ) 
the three Pauli matrices. Moreover, < X < 00, 
< 9 < it/2, < <p 1 < 2tt, and < ip 2 < 2tt. TRS 
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<7, 



s 1 



cr 2 
cr 9 



(2.2) 



The matrix S is the most general 4x4 unitary matrix 
that describes a quantum tunneling process between two 
Kramers doublets that preserves TRS but can break SRS. 
When r = 1 and t = 0, S is reduced to the unit ma- 
trix, and there is no tunneling between one Kramers pair 
(tp 1 + ip 2 ) and the other pair (tp 3 + ip A ). The tunnel- 
ing with (without) a spin flip occurs with the probability 
t 2 sin 2 9 (t 2 cos 9) . Although S is parametrized by four 
real parameters, only X and 9 matter as ip 1 and ip 2 can be 
absorbed in the overall phase windings that the Kramers 
doublets acquire when traversing along links connecting 
nodes. 

Disorder is introduced in the network model by assum- 
ing that the phases of Kramers doublets on the links and 
the phases 9 on the nodes are independently and iden- 
tically distributed. The distribution of the link phase of 
a Kramers doublet is uniform over the interval [0, 27r[. 
The distribution of 9 is sin 29 over the interval [0,7r/2]. 
We are left with one parameter X in the network model 
that controls the scattering amplitude at every nodes. 
The parameter 9 in the network model plays the same 
role as spin-orbit interactions (of Rashba type) in a ran- 
dom tight-binding Hamiltonian belonging to the two- 
dimensional symplectic symmetry class, whereas the pa- 
rameter X plays the role of the Fermi energy. 

In order to distinguish the topologically-trivial insulat- 
ing phase from the QSH insulating phase, we shall com- 
pare the results that we obtained from the network model 
against the ones that we obtained from a two-dimensional 
tight-binding model introduced in Ref. ||, the so-called 
SU(2) model. The SU(2) model is a microscopic random 
tight-binding Hamiltonian with on-site randomness (box 
distributed with the width W) and with random hop- 
ping amplitudes such that the spin-dependent hopping 
amplitudes between any pair of nearest-neighbor sites 
of a square lattice are distributed so as to generate the 
SU(2) invariant Haar measure. It captures the transition 
between a metallic and a topologically-trivial insulating 
state in the ordinary two-dimensional symplectic univer- 
sality class. In the SU(2) model, the Fermi energy plays 
the role of the network parameter X for a fixed and not 
too strong W . 



III. PHASE DIAGRAM AND BOUNDARY 
CONDITIONS 



On symmetry grounds, we expect that it is possible to 
drive the network model through two successive Ander- 
son transitions by tuning X. Indeed, this was shown to be 
the case in Ref. 12. For any X bounded by the two critical 
values X s < X l [X s = 0.047±0.001, X l = 0.971±0.001; 
the subscript s (I) stands for small (large)] the network 
model is in a metallic state, while for X < X„ or for 
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FIG. 2: (Color online) Phase diagram for the network model 
after Ref. |l|. 



X l < X it is in an insulating state. It was shown in 
Ref. [12] that the localization length, a bulk property of 
any one of the two insulating phases, diverges with the 
exponent v ps 2.7 of thc-ordinary two-dimensional sym- 
plectic universality classcl upon approaching the mobility 
edge (either X s or X t ). The resulting phase diagram 
is shown in Fig. |[ The conclusions from Ref. [l^ were 
reached by studying numerically the network model in a 
strip geometry with periodic boundary conditions (PBC) 
imposed in the transverse direction. 

Here, we shall depart from Ref. [l| by imposing reflect- 
ing boundary conditions (RBC) in the transverse direc- 
tion in the strip geometry such as in Fig. or Fig. 0(a) for 
example. Reflecting boundaries for the network model 
that respect the TRS arc defined by choosing a set of 
nodes and amputating two pairs of links (one labeled by 
an odd integer and the other by an even integer in the 
convention of Fig. 0) attached to any of these nodes in a 
consistent fashion, i.e., no node can be the end point of 
an odd number of pairs of links. A boundary node is de- 
picted by a colored semi-circle, and its scattering matrix 
S is always represented by a unit 2x2 matrix, regard- 
less of its index S or S'. Three examples of reflecting 
boundaries defining a strip geometry are shown in Fig. 0. 
From the fact that the scattering process in Fig. 0(b) is 
obtained by a clockwise rotation by tt/2 of the scattering 
process in Fig. 0(a) follows the important property that 
a horizontal boundary passing through nodes of type S 
is equivalent to a vertical boundary passing through the 
nodes of type S and vice versa. 

In the following discussions wc shall focus on the crit- 
ical point X — X x and its neighboring insulating phase. 
Imposing RBC, we shall sec that the insulating phase at 

X > X[ in the phase diagram from Fig. can acquire a 
topological attribute in that it supports a single Kramers 
doublet that is localized in the direction transverse to the 
boundary but is delocalized along the boundary, i.e., a 
single Kramers degenerate pair of edge states. 

The refined phase diagram Fig. is depicted for X > 

X I according to the presence or absence of edge states 
in the insulating phase at X > X t . Identifying helical 
edge states can be done in a pictorial way in the perfect 
reflection limit X — > 00, provided we assume continuity 
of the insulating phase. In the limit X — > 00, the scatter- 
ing at the nodes S and S' reduces to the glancing events 
depicted in Fig. 0. As a corollary Fig. ]| simplifies to 
Fig. 0. We thus see two Kramers doublets propagating 
along the upper and lower boundaries in Fig. 0(a), none 
111 Fig. 0(b), and one in Fig. 0(c). The refined phase-dia- 
grams from Fig. then follow assuming continuity.E-3 In 
the following sections we go beyond the pictorial argu- 





FIG. 3: (Color online) Quasi-one-dimensional network model 
with reflecting boundaries, (a) The upper and lower bound- 
aries pass through the nodes S'. (b) The upper and lower 
boundaries pass through the nodes S. (c) The upper and 
lower boundaries pass through the nodes S and S', respec- 
tively. Transverse periodic boundary conditions are only pos- 
sible in geometries (a) or (b) for which the transverse width 
is here M — 2. 



ment from Fig. and study critical properties at the An- 
derson metal-insulator transition in a quantitative fash- 
ion with an emphasis on the boundary and corner multi- 
fractal scaling exponents at criticality. 

It is important to note here that we may consider the 
insulating state in Fig. 0(a) as a Z 2 topological insulating 
state because of the presence of a helical edge state along 
the boundaries. Similarly, we may identify the insulating 
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FIG. 4: (a) The network model on a square with L nodes 
of type S and L 2 nodes of type S' (a node on a boundary is 
counted as 1/2) is wrapped around a cylinder by imposing 
reflecting boundary conditions (RBC) along the x direction 
and periodic boundary conditions (PBC) along the y direc- 
tion, (b) It is wrapped around a torus by imposing PBC in 
the x and y directions. 
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FIG. 5: Phase diagram of the network model of the quasi- 
one- dimensional geometries in Fig. |^(a), |^(b), and [](c), re- 
spectively, at X > X l . The choice of the boundary decides 
whether the insulating phase at X > X t has helical edge 
states. 



state in Fig. ||(b) with a Z 2 topologically trivial, ordinary 
insulator. These observations indicate that, in the net- 
work model formulation, the Z 2 topological nature of an 
insulating state can be determined by appropriate choice 
of reflecting boundary nodes, even without changing the 
control parameter X. This is very similar to the situa- 
tion in the Chalker-Coddington network modclO for an 
integer quantum Hall plateau transition, in which an in- 
sulating state may or may not have an edge state for a 
given value of quantum tunneling parameter at nodes, 
depending on the location of reflecting boundary nodes; 
an insulating state with an edge state can be regarded 
as an integer quantum Hall state. We also note that we 
may attribute this dependence of the topological nature 
of an insulator on the reflecting boundary conditions in 
network models to our freedom to give any topological 
character to the vacuum outside the insulator. 
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FIG. 6: (Color online) Scattering in the reflective limit X 
oo for node S (a) and S' (b). 



IV. 



NORMALIZED LOCALIZATION LENGTH 



Supporting evidences for the phase diagram in 
Fig. ||(a) can be extracted from the dependence on M 
of 



A (t) {X,M)=t (l) {X,M)/M 



(4.1) 



in the geometry of Fig. ||(a). In Eq. (4.1), the ith normal- 
ized localization length hy'(X,M) is the ratio between 
the ith localization length ^'(X, M), which is-given by 
the inverse of the value of the ith smallest pairEJ of Lya- 
punov exponent of the transfer matrix, and the width M. 
in either one of the geometries of Figs. ||(a) or |](b).EZl 
The transfer matrix for the strip geometries can be con- 
structed in a similpj manner as in the cylinder geometry 
(transverse PBC).tL3 

Figure || shows the dependence on X of A' 1 ' (X, M) 
with i = l and i = 2 for the values of M = 4, 8, 
and 16, respectively, in the geometry of Fig. ||(a). Also 
shown in Fig. |^ is the dependence of As l '(X,M) for 
transverse PBC. The existence of the critical points X s 
and X[ is signaled by the independence of A.W(X t , M) 
on M. Imposing transverse PBC reduces the finite-size 
corrections and allows a more accurate measurement of 
X s and X l , the dashed vertical lines in Fig. ||Ji3 When 
X s < X < Xi, AW(Jf,M) increases with M for both 
transverse RBC and transverse PBC, as is expected from 
a two-dimensional metallic state. This trend is reversed 
in the insulating phases for transverse PBC, as is ex- 
pected from a two-dimensional insulating state. How- 
ever, when X l < X and transverse RBC are chosen, 
A*4'(X, M) remains an increasing function of M whereas 
A( 2 '(X, M) becomes a decreasing function of M (the in- 
set shows this more clearly with its logarithmic vertical 
scale). This opposite dependence on M of A^(X,M) 
and A^ 2 ) (X, M) is the signal that each of the upper and 
lower boundaries supports a single Kraniers_doublet that 
would be extended along the boundary,t§t£l were it not 
for the existence of a finite tunneling amplitude (i) that 
couples the two Kramers doublets (edge states) residing 
near the two opposite boundaries and (ii) that is expo- 
nentially small in M when M is much longer than the 
mean free path. The dependence on M of K^{X,M) 
with transverse RBC when X < X s is the one expected 
from an ordinary insulating state. All together, these 
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FIG. 7: (Color online) Large X limit of the quasi-one- 
dimensional network model with the reflecting boundaries 
from Fig. H. 



observations point to the refined phase diagram shown 
in Fig. ||(a), i.e., the insulating phase is unconventional 
when X l < X due to the presence of a single Kramers 
doublet of edge states per boundary that becomes delo- 
calized along the boundaries in the limit M — > oo. 

The M- independent, normalized localization length at 
the Anderson transition is nothing but the normalized 
correlation length A c at a generic continuous phase tran- 
sition in the theory of critical phenomena. This quan- 
tity is known to depend on the choice of the transverse 
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FIG. 8: (Color online) The dependence on X of A (1) (X, M) 
in the geometry of Fig. ^|(a) with transverse RBC is shown 
for M = 4 (0), 8 (■), and 16 (•). The dependence on X of 
A (2) (X, M) with transverse RBC is shown for M = 4 (x), 8 
(+), and 16 (*). The dependence on X of A W (X,M) with 
transverse PBC is shown for M = 4 (blue dashed curve), 8 
(green dotted curve), and 16 (red solid curve). The vertical 
dashed lines identify the critical points X B and Xi deduced in 
Ref. [H^ when transverse PBC are imposed. The inset displays 
the same data points with a logarithmic vertical scale. 



boundary conditions .0 First, the value 

A (pbc) (tT c ) = 1.844 ±0.002 



(4.2) 



of the normalized localization length at the Anderson 
transition was calculated in Ref. ^ for the SU(2) model 
with transverse PBC. It should here be compared with 
the values 



A (pbc)( X j = L81±aolj 
A (pbc) (A ; ) = 1.82 ±0.01, 



(4.3) 



calculated in Ref. 12|. Second, we have extended the cal- 
culation in Ref. |2 for the SU(2) model in which trans- 
verse PBC were used to the case in which the vanishing 
of the wave functions along the transverse boundaries is 
imposed, a situation that we shall refer to as fixed bound- 
ary conditions (FBC), and found [see Fig. 0(c)] 



A (fbc) (lT c ) = 1.50 ±0.03. 
This value should be compared with 

A (rbc) (Aj = 1.49 ±0.02 



(4.4) 



(4.5) 



for the network model with transverse RBC [where we 
can also deduce the right-hand side from Fig. |9|(c) for 
which it is the network model at the critical point X, in 
the geometry of Fig. ||(b) that is investigated; Fig. |](a) 
shows insulating behavior of A^ for X > A ; ] . Evidently, 
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FIG. 9: (Color online) (a) Dependence on X of the nor- 
malized localization length A (1) (X, M) = A(X,M) for the 
network model in the quasi-one-dimensional geometry of 
Fig. ||(b) with transverse RBC and M = 4, 8, 16, 32, 64. (b) 
Dependence on the dimensionless disorder strength W of 
the normalized localization length A W (W,M) = A(W,M) 
for the SU(2) model in the quasi-one-dimensional geome- 
try with transverse fixed boundary conditions (FBC) and 
M — 4, 8, 16, 32, 64. The solid curves are computed from the 
finite-size scaling functions that follow, (c) Finite-size scaling 
analysis for the data shown in (a-b): The solid blue curves 
represent the scaling functions for the network model with 
RBC on node S in the vicinity of the critical point X t where 
we find A (rbc) (A) ) = 1.49 ±0.02 and v = 2.88 ±0.04. The red 
dashed curves represent the scaling functions for the SU(2) 
model with FBC in the vicinity of the critical point W c where 
we find A (rbc) (W c ) = 1.50 ± 0.03 and v = 2.85 ± 0.06. The 
scaling functions are obtained from finite-size scaling analy- 
sis incorporatjipg corrections from a leading irrelevant scal- 
ing variable.li3'EJ The normalized localization length A' is ob- 
tained from A by subtracting these corrections, as defined in 
Eqs. (3.7) and (3.8) from Ref. ||. The distance to the critical 
point \X — Xi is rescaled by a factor c ~ 1.7 in the scaling 
function for the network model with RBC so that it coincides 
with that for the SU(2) model with FBC. 



the values (4.4) and (|4.5|) agree within their error bars 



(in support of our identification, made at the end of the 
last section, of the insulating phase without a helical edge 
state with a Z 2 topologically trivi al in sulator) but clearly 
differ from the values (4.2) and (4.3). Finally, and more 
importantly, there is a rather large asymmetry 



A( rbc )(X ; )/A< rbc )(A: s ) pa 4.8 



(4.6) 



for the network model with transverse RBC in the geom- 
ctry of Fig. |(a). 

Finite-size scaling in Fig. ^|(c) clearly shows that the 
scaling function obtained for the SU(2) model with trans- 



verse FBC is identical to that for the network model at 
X = X l in the geometry of Fig. ||(b). These scaling 
functions are therefore universal property of the critical 
point between a metallic phase and an ordinary insula- 
tor. From this finite-scaling analysis we also obtained the 
critical exponent v of the diverging localization length. 
At the critical point X x of the network model we found 
v = 2.88 ± 0.04. Again, this should be compared with 
the critical exponent v of the SU(2) model with trans- 
verse FBC for which we find v = 2.85 ± 0.06. Both 
exponents agree with each other within their error bars; 
they are also consistent with the-exponent 2.7 < v < 2.8 
obtained with transverse PBC.cni3 This implies that the 
exponent v is a bulk property independent of boundary 
conditions (whereas the scaling functions are dependent 
on the transverse boundary conditions). 

The network and SU(2) models in a quasi-one- 
dimensional cylinder geometry at criticality are indis- 
tinguishable as measured by the normalized localization 
length (4.2) and (fD3), respectively. The same is true 
in a strip geometry, provided the insulating side of the 
transition in the network model is topologically tr ivia l 
(i.e., has no helical edge states) according to Eqs. (4.4) 
and (4.5). However, the strong asymmetry ((Oj) hints 



at the possibility that some boundary critical exponents 
might be sensitive to the choice of transverse boundary 
conditions that dictates the presence or absence of a de- 
localized Kramers doublet of edge states. 



V. BULK AND BOUNDARY CRITICALITY 
A. Typical spatial profile of wave functions 

To gain more insights on the criticality of the net- 
work model, we study numerically the critical normal- 
ized wave functions \I r when X = X l in the geometries 
of Figs. |i|(a) or f|(b). To this end, the support of nor- 
malized wave functions for the network model is defined 
on the midpoints (x, y) of every bond joining nearest- 
neighbor nodes, i.e., the normalized wave function \1/ can 
be viewed as a complex-valued vector whose 8L 2 com- 
ponents ip a (x,y), with the spin-1/2 quantum number la- 
beled by a =T,i, correspond to 4L 2 freely propagating 
Kramers pairs of plane waves. The dynamics of such a 
wave function is governed by a unitary evolution operator 
U built o-Hjt of all scattering matrices at the nodes of the 
network.E2l For each realization of the disorder, wc nu- 
merically diagonalize U and retain one normalized wave 
function whose eigenvalue is closest to 1. The number 
of disorder realizations is 5 x 10 4 for each system size 
L (L ranging from 20 to 80). By tuning the parameter 
X = X l at each node, this normalized wave function is, 
with a slight abuse of language, called critical. 

First, wc define 
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metal and an ordinary insulator. Remarkably, Fig. |To|(b) 
shows that the dependence on x of (\xi\^\ 2 } x L is non- 
monotonic when approaching a boundary from the cen- 
ter of the bulk x — L and that its absolute maximum is 
reached at the boundaries x = 1, 2L instead of the local 
maximum at x = L. 



B. Scaling exponents A q 

This striking difference between Fig. |R](a) and 
Fig. |l^(b) suggests that the statistics of the critical nor- 
malized wave functions is sensitive to the topological na- 
ture (i.e., the presence or absence of helical edge states) 
of the insulating state. A powerful tool to study the 
statistics of critical normalized wave functions is the mul- 
tifractal scaling analysis Multifractality is encoded by 
a set of scaling exponents (anomalous dimensions) A q K ^ 
defined by the scaling laws 
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FIG. 10: (Color online) Dependence on the position x 



Eq. (|y) for the values of L = 20, 25, 30, 35, 40, 50, 60, 80 from 
top to bottom with (a) S boundaries or with (b) S' boundaries. 
The blue dashed curves in (a) represent the normalized crit- 
ical wave functions (lnl*! 2 )^ L for L = 36,48,60,72,96, 120 
obtained from the SU(2) model defined in Ref. ^ using the 
cylinder geometry of Fig. M(a). 



for L > 1, where ^ is the normalized spinor obtained 
from the network model in the geometries ^(a) or [|(b) 
along the cylinder axis of Fig. h|a) of (ln]*! 2 )^ defined by a t criticality X = A,, and q is any real number. The 

index (k) distinguishes bulk, boundary, and corner ex- 
ponents through the choice of the location of r = (x,y) 
relatbic_tp^lJicJpoundaries and the nature of the bound- 
aries. BOoo'E£l In this paper, we distinguish six cases, 
(k) = (2), (1,Z 2 ), (1,0), (0,Z 2 ), (0,O), and (0,Z 2 |O). 
Here the first entry in (k) refers to the dimensionality of 
the support of wave functions. The second entry in (k), 
Z 2 or O ("ordinary"), refers to the presence or absence, 
on the insulating side of the critical point, of a Kramers 
doublet (helical) edge state al ong a boundary or at a cor- 
ner where the measurement (5^) is made (see below for 
more deta iled definitions). 

In Eq. ( |5.2] ) the proportionality constant may have a 
dependence on L that is much milder than the power law 

it multiplies. This means that the scaling exponents A q K ^ 
are obtained from calculating numerically 



where X is tuned to the critical point X l for the network 
model. The ovcrline denotes disorder averaging over the 
nodes through the independently distributed angles 9 as 
well as over the links through the independently and uni- 
formly distributed random phases of Kramers doublets. 

Figure [h] shows the x dependence of (lnl^l 2 ^ . for 
different values of L at the critical point X l in the cylin- 
der geometry of Fig. ||(a). The two figures correspond 
to different RBC in the x direction: (a) the boundaries 
passing through the S nodes as in Fig. ||(b), and (b) the 
boundaries passing through the S' nodes as in Fig. |3|(a). 
The boundary bonds are located at x = 1 and x = 2L. At 
the critical point between the metallic and the insulating 
state without helical edge state [Fig. |T^(a)] , the depen- 
dence of (in | 'I' | 2 )^ L on x is symmetric about x = L 
where it reaches a maximum, and reaches a minimum 
close to the boundaries. For comparison, we show in 
Fig. ^(a) (In l*! 2 )^ L of the SU(2) model, in which the 

on-site disorder strength W i s tun ed to the critical point 
W c , and the overlinc in Eq. ( |5.lD denotes disorder aver- 
aging over the on-site energies and the Haar measure of 
the hopping matrix elements. We find no qualitative dif- 
ference in the dependence on x and Lof(ln|\E , | 2 ) ;cL be- 
tween the two models when the critical point separates a 



D(K) ^ gln|vl/(r)| 2 -ln|vl/(r) 
q hi L 



\2q 



(5.3) 



first. Here, the double overline means that in addition to 
the disorder averaging a spatial average over the relevant 
sites r (AL 2 sites for the torus geometry, 2L sites for the 
boundary) is also taken to improve the statistics. This is 

followed by a linear fit of (L) as a function of 1/ In L. 
Only system sizes larger than L = 35 are kept for the 
linear fit. At last, Aq is obtained as the intercept of 
the linear fit with the vertical axis at 1/lnL = 0. For 
the case of boundary scaling exponents, we also replaced 
^(r)! 2 with r a bond joining a node on the boundary by 
averaging it over all bonds from the elementary plaquctte 
to which it belongs. This coarse graining is necessary to 
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FIG. 11: (Color online) The dependence on q of A^' /[q{l — 
q)] (•) when X = X l for (a) the bulk: (k) = (2), (b) the 
S boundary: («) = (1,0), and (c) the S' boundary: (/t) = 
(1,Z 2 ). The solid line s rep resent a^ K \q) — 2 evaluated at 
q — according to Eq. (5.4a) with its error bars indicated by 
dashed lines. The dependence on q of Ag K '/[g(l — g)] for the 
SU(2) model in the bulk and on the boundary are shown with 
o in (a) and (b), respectively. Blue open circles o in (c) show 
the q dependence of A^ 2 '/[q(l — q)]. 



overcome the small oscillations that are visible in the x 
dependence plotted in Fig. [To] close to the boundaries. 

Figure fir] shows the q dependence of the scaling ex- 
ponents Alp for the network and the SU(2) model in 
different geometries. 

The bulk scaling exponent s, fo r which case (k) = (2), 
characterize the scaling law ( |5.2| ) in the torus geometry 
from Fig. ^(b) (which has no boundary). This is also the 
case in the cylinder geometry of Fig. |](a), if the distance 
from the site r to the boundaries is of order L. The 
bulk exponents shown in Fig. |ll|(a) are obtained with 
the torus geometry in order to be free from boundary 
effects and to maximize the number of sampling points r. 
Figure |ll](a) shows that the bulk scaling exponents A q 2 ^ 
of the network model are identical to those obtained for 
the SU(2) model] within their error bars. 

The boundary scaling exponents are obtained by 
putting r on the boundaries (with the coarse grain- 
ing mentioned above) in the cylinder geometry from 



Fig. |J(a). There is an important caveat however. As we 
saw in the preceding sections, the network model sup- 
ports two types of (straight) boundaries, while the SU(2) 
model has only one. Indeed, for the network model at 
X = X t studied numerically in Fig. [ll]we have two cases: 

(i) (k) = (1,0), when the boundary on which r is lo- 
cated passes through the nodes of type S. To reduce 
statistical errors, we take the cylinder geometry of 
Fig. ^(a) which follows from imposing longitudinal 
PBC in Fig. ||(b). No helical edge state exists on 
the boundaries in the insulating phase at X > X l . 

(ii) (k) = (1,Z 2 ), when the boundary on which r is 
located passes through the nodes of type S'. To 
reduce statistical errors, we take the cylinder ge- 
ometry of Fig. |](a) which follows from imposing 
longitudinal PBC in Fig. [|(a). There exists a heli- 
cal edge mode on each boundary in the insulating 
phase at X > X l . 

Figure |ll|(b) shows that the boundary scaling exponents 
Aq 1 ' ' in the network model agree, within their error 
bars, with the boundary scaling exponents obtained in 
Ref. H from the SU(2) model in the cylinder geometry 
of Fig. [|(a). Remarkably, Fig. |ll|(c) shows a second 
set of boundary scaling exponents Aq 1 ' 1 " 2 ^ in the net- 
work model, which is markedly different from the first 
set of boundary scaling exponents Ag 1,0 ' in Fig. |TT|(b) . 
To sum up, the three distinct sets of scaling exponents 



A 



(2) A (l,0) 



, and Aq 2 are identified from Figs. pd|(a), 
|ll|(b), and |n](c), respectively. 

The emerging picture is that the network and the 
SU(2) models at criticality share common bulk scaling 

(2) 

exponents A q and a set of boundary scaling exponents 
Aq at an ordinary boundary that docs not support 
a helical edge mode in the neighboring insulating phase. 
Only the network model has one more set of boundary 
scaling exponents Aq 1 ' 1 " 2 at a boundary that has a helical 
edge mode in the insulating phase. 



C. /(a) spectra 

There is yet another set of scaling exponents f^ K \a^) 
that we have calculated for the SU(2) and network mod- 
els. These exponents are calculated numerically from the 
scaling ansatz 



*(r)| 2g ln|*(r)|" 



Mr) 



2g 



~a q K) hxL 



(5.4a) 



and 



ln|*(r) 



2q 



/(")(«(*)) - a g K) q - d M ] lni, (5.4b) 



where = 2 and 1 for the bulk and boundary expo- 
nents, respectively. It can be shown that the Legcndrc 
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FIG. 12: (Color online) (a) Dependence on aS K ' of the multi- 
fractal spectra in the bulk (k) = (2) (a), at the S bound- 
aries (k) = (1, O) (0), and at the S' boundaries («) = (1,Z 2 ) 
(•), at X — X[ . Dependencies on q of a q and f q are shown in 
(b) and (c), respectively, with the same symbols as in (a). The 
bulk and boundary multifractal spectra for the SU(2) model 
defined in Ref. ^ are shown by solid and dashed curves, re- 
spectively. 



transforrJ@0ffl 

/ (K) (a (K) ) = (a (K) - 2) 9 - A[ k} + d {K) , (5.5a) 
where q is a function of a^ K ' obtained from inverting 

a^(q)-2 = dA^/dq, (5.5b) 



relates the scaling exponents (5.2) and (fo-4). The number 
y(« ; )( Q! ( K )) i s the fractal (i.e., HausdorfEJ) dimension of 

the set of points r such that |\T/(r)| 2 scales as L~ aM . 

The dependence on of for the critical network 
model at X = X l in the torus, (k) = (2), in the cylinder 
geometry with S boundaries [Fig. 0(b)], (k) = (1, O), and 
in the cylinder geometry with S' boundaries [Fig. p|(a)], 
(k) = (1, Z 2 ), are shown in Fig. |l2|(a), which are obtained 
by combining a q in Fig. |T^(b) and fq^ in Fig. |l^(c). 
We see that there are three distinct multifractal spectra 
at the critical point X l : for the bulk, /t 1 ' ) for the S 
cylindrical geometry, and f^ 1 ^) for the S' cylindrical ge- 
ometry. They are compared with the multifractal spectra 



for the SU(2) model in the torus and cylinder geometry 
obtained in Ref. [j| Within their error bars, they agree 
with /( 2 ) and f^ l '°\ respectiv ely. The exponent a q at 
q = is closely related to Eq. (5.1) with x chosen in the 
bulk or on a boundary. For the network model we find 
the values 



a 



(2) 




The inequality oiq 1,2 



2.174 ±0.001, 
= 2.420 ±0.005, 
= 2.086 ±0.015. 



(5.6a) 
(5.6b) 
(5.6c) 



< a 



(2) 




(l.O) 
*0 



the x dependence seen at the boundaries in Fig. |10|, as 
wave functions near S' boundary in the S' cylindrical ge- 
ometry are expected to be more extended because of the 
existence of edge modes in the insulating side. We note 
that the error bars are an order of magnitude larger for 
ay because of the presence of larger finit e-size cor- 
rections. Within their error bars, the values ( [5. 6a ) and 
( |5.6b| ) agree with the-pOnes for the SU(2) model in the 
bulk and boundaries,™ respectively. 

It is worth mentioning that the multifractal analysis 
performed here involves extracting scaling exponents af- 
ter performing the disorder averaging. Extracting scal- 
ing exponents before performing the disorder averaging 
yields typical scaling exponents. Typical scaling expo- 
nents need not be id en t ical .with average scaling expo- 
nents calculated hereB3'E2(E2l The average scaling expo- 
nents Aq are expected to differ from the typical-ones 
for any values of q such that f^(a^ K ^ ) is negative;^ see 
Fig. pL2|(c) . For this range of q rare events dominate the 
calculation of A q as is evidenced by the larger error bars 
in Fig. This explains the systematic deviations from 
the mirror symmetry^ 2 ] about q = 1/2 of A q K ^ /[q(q — 1)] 
for large q in Fig. O. 

The rational for studying these average scaling expo- 
nents is that they are expected to be the scaling di- 
mensions of some primary operators representing mo- 
ments of wave function amplitudes^ an underlying two- 
dimensional conformal field theory£3 Knowing them con- 
strains the possible field theories that can encode critical 
properties of an Anderson transition. 



VI. CORNER MULTIFRACTALITY 

So far we have always considered geometries of the 
network model with boundaries of the same type, as in 
Figs. ||(a) or ||(b). We are now going to investigate the 
case of mixed boundaries. 

We begin with the network model in the quasi-one- 
dimensional geometry (quantum wire geometry) from 
Fig. ||(c) for which we show that the Landauer con- 
ductance obeys a statistical distribution that differs 
from that of an ordinary symplectic quantum wire, i.e., 
a symplectic quantum wire in a geometry compatible 
with transverse PBC. A physical realization of this case 
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FIG. 13: (Color online) (a) A semi-infinite geometry with two 
point contacts (green dot-dashed curves) attached at the two 
interfaces between different types of boundaries. Each line 
or curve represents a Kramers doublet, (b) Closed network 
with mixed boundaries. The thick wavy and solid lines on 
the edges represent two different types of boundaries, with 
and without a helical edge mode at X > X L , respectively. 



A. Quasi-one-dimensional wire with mixed 
transverse open boundary conditions 

In the quasi-one-dimensional wire depicted in Fig. |^(c), 
whereby the length of the disordered region N is larger 
than the mean free path t 1, one transmission eigen- 
value is always unity for any disorder configuration while 
the other 2M eigenvalues are exponentially small with 
N. This follows from applying the results from Refs. |34| 
and |l9[ The same is true of the thick quantum wire limit 
N, M — > oo with the ratio Ml/N held fixed, as it is now 
the analysis from Rcf. [l8| that can be borrowed. 

The persistence of the unit eigenvalue of the trans- 
mission matrix is unique to the network model with 
mixed boundary conditions. This remarkable property 
can be understood physically as follows. Because of the 
mixed boundaries, one and only one edge supports a sin- 
gle Kramers doublet of edge states. This Kramers dou- 
blet of edge states is free to propagate along the edge 
and thus provide a single perfectly conducting channel. 
This physics is indeed realized with helical edge modes 
of a two-dimensional quantum spin Hall insulator and 
a metallic carbon nanotube in the absence of interval- 
ley scattering. Mathematically, the transfer matrix A4 
for the geometty of Fig. ||(c) is a member of the group 
SO*(4M + 2).E3 The eigenvalues of ln{MM^) are sym- 
metrically distributed around and denoted by ±2x i 
(0 < < x 2 < ••■ < x 2M+i)- I n turn, the trans- 
mission eigenvalues T i are written as T i =1/ cosh 2 x i 
(i = 1,2, ... , 2M + 1). Because of the Kramers degen- 
eracy, the eigenvalues of h^A^A'T) are two-fold degener- 
ate. This leads to the conclusion that there must be a 
null eigenvalue x± = 0, thereby T x = 1. 

B. Two point contacts on a boundary of mixed 
type 



is given by metallic carbon nanotubes with spatially 
smooth disorder potential. As shown by Ando and 
SuzuuraJ^J in the absence of intervalley scattering by dis- 
order, conduction of electronic states near a Fermi point 
is described by transfer matrices in the symplectic class 
with an odd number of conduction channels. A mathe- 
matical treatment—of such conduction process was initu 
ated by Zirnbaucnia and further developed by Takane.Ea 

We then turn our attention to a semi-infinite geometry 
with two point contacts to ideal reservoirs attached at the 
interface between the boundary of type S and of type S' 
as shown in Fig. [l^(a). After tuning the network model 
to criticality, we shall see that this amounts to studying 
corner multifractal scaling exponents whose values differ 
from those taken by the corner multifractal scaling at the 
corner along a homogeneous segment of a boundary, see 
Fig. [11(b). 

Common to both examples of mixed boundaries is the 
existence of a single perfectly conducting channel. 



Next, we consider the case of a semi-infinite network 
model with a boundary that is mixed, whereby it is nec- 
essary to attach two point contacts (leads) at the two in- 
terfaces between the boundary nodes of type S and S' in 
order to maintain TRS. This setup is shown in Fig. p~3|(a) . 
The scattering matrix that relates incoming to outgoing 
waves from the leads is then a 2 x 2 matrix. The con- 
straint that it belongs to the symplectic symmetry class 
forces this matrix to be proportional to the unit 2x2 
matrix up to an overall (random) phase. Hence, the two- 
point conductance is always unity however far separated 
the two point contacts are. An incoming Kramers dou- 
blet is transmitted with probability one through the dis- 
ordered region in Fig. |l|(b). This result may be related 
to the presence of a perfect transmission channel in the 
geometry of Fig. [|(c) by thinking of a conformal map- 
ping that transforms a half plane with mixed boundaries 
to an infinite strip with mixed boundaries. 

Before proceeding with the case at hand, we recall that 
it is expected on general grounds that the moments of the 
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two-point conductance in a network model at criticality, 
when the point contacts are far apart, decay as power 
laws with scaling exponents proportional to the scaling 
exponents A^.EJ Consequently, after tuning the semi- 
infinite network model to criticality, the zero-dimensional 
boundary scaling exponents for the moments of ^(r)) 2 
with r at an interface of different types of boundaries 



fo,z 3 |o) 







(6.1) 



emerge as a signature of the nontrivial topological nature 
of the insulating side at the Anderson transition. We con- 
jecture that, if a description of the critical point e xist s 
in terms of a conformal field theory, the exponents (6.1) 
might then be obtained from the correlation functions be- 
tween the moments of the local operator encoding the lo- 
cal density of states and additional-insertions of a bound- 
ary condition changing operator .L3 

The zero-dimensional multifractal scaling expo- 



nents (6.1) are different from the corner multifractal scal- 
ing exponents which characterize scaling of the moments 
of wave functions at a corner with the boundary of a given 
type, such as the upper left corner and the lower right 
corner in Fig. |l3|(b). These corner multifractal exponents 
read 



VII. CONCLUSIONS 



We have shown that the boundary multifractal spectra 
for the critical normalized wave functions are sensitive to 
the choice of boundary conditions at the Anderson tran- 
sition in the two-dimensional symplectic class. This is 
the first example where boundary multifractal exponents 
are calculated under different boundary conditions in the 
problem of Anderson localization-delocalization transi- 
tion. It would be interesting to look for other examples 
of disorder-induced continuous phase transitions at which 
boundary critical properties depend on boundary condi- 
tions but bulk critical properties do not. 

We conjecture that the two-dimensional conformal 
theory describing the Anderson transition in the two- 
dimensional symplectic universality class, if it exists, 
should be compatible with two distinct conformally in- 
variant boundary conditions on the boundary of a half- 
plane, thereby yielding two distinct sets of boundary mul- 



The recent obser- 



tifractal exponents A^ 1 ' ^ and A, 
vatioja of the QSH effect in HgTe/(Hg,Cd)Te quantum 
wellsO suggests that it might be possible to probe exper- 
imentally the Kramers degenerate edge states at critical- 
ity and the corresponding boundary multifractal spectra. 



a(o,o) 
q 



2A(!.o) 



(6.2) 



A 



(0,Z, 



2A 



(i, 



(6.3) 



depending on the type of the boundary and the criti- 
cal point to which the network model has been tuned. 



Here, Eqs. (3.2) and (6.3) follow from a general relation 



based on two-dimensional conformal mappings that re- 
lates the scaling exponents ( |5.2[ ) with r at the corner 
between two straight boundaries m eeti ng with the an- 
gle 7r/2 and the scaling exponents (5.2) with r along a 
straight boundaryE 
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